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Abstract 

Existing  analysis  of  the  carrier  sense  multiple  access  mode 
(CSMA)  has  led  to  the  determination  of  the  average  channel  performance 
in  terms  of  average  throughput  and  average  packet  delay.  This  was 
achieved  by  formulating  a  semi-markovian  model  for  CSMA  channels  with 
a  finite  population  of  users  [1,2].  In  this  paper,  it  is  shown  that, 
using  the  same  model,  it  is  possible  to  derive  the  actual  packet 
delay  distribution.  The  analysis  is  similar  in  nature  to  that  pro¬ 
vided  for  slotted  ALOHA  channels  in  [3].  These  results  are  useful 
in  the  design  of  systems  intended  fcr  real  time  applications  such 
as  digitized  speech,  and  in  the  analysis  of  multihop  packet  radio 
systems . 
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I .  INTRODUCTION 


Carrier  sense  multiple  access  (CSMA)  is  an  efficient  random 
access  method  for  multiplexing  a  population  of  users  communicating 
over  a  shared  packet-switched  channel.  Briefly,  CSMA  reduces  the  level 
of  interference  (caused  by  overlapping  packets)  in  the  random  multi¬ 
access  environment  by  allowing  terminals  to  sense  the  carrier  due 
to  other  users'  transmissions;  based  on  this  channel  state  infor¬ 
mation  (busy  or  idle) ,  the  terminal  takes  an  action  prescribed  by 
the  particular  CSMA  protocol  in  use.  In  particular,  terminals  never 
transmit  when  they  sense  that  the  channel  is  busy  [A]. 

To  analyze  the  performance  of  CSMA,  a  semi-markovian  model 
has  been  formulated  for  channels  with  finite  populations  of  inter¬ 
active  users  [1,2].  However  only  average  stationary  performance 
was  derived  in  terms  of  average  throughput  and  average  packet  delay. 
Unfortunately,  average  performance  may  not  be  adequate  when  designing 
systems  intended  for  real  time  applications  such  as  digitized  speech, 
and  the  analysis  has  to  be  extended  so  as  to  include  delay  distri¬ 
butions.  Also,  when  analyzing  multihop  systems,  it  is  important 
to  be  able  to  characterize  the  departure  process  from  a  collection 
of  nodes,  as  this  corresponds  to  the  arrival  process  to  other  nodes 
in  the  multihop  environment. 

In  this  paper,  we  extend  this  analysis  and  show  that, 
using  the  same  model,  one  can  derive  the  actual  distribution  of 
packet  delay  and  interdeparture  time.  In  section  II,  we  describe 
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^  - 

the  model.  In  section  III,  we  review  the  derivation  of  the  average 
channel  performance  as  presented  in  [1,2].  In  section  IV  we  address 
the.  interdeparture  time  distribution  and  its  moments.  In  section  V, 
we  treat  packet  delay. 
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II.  THE  MODEL 


Although  the  operation  of  CSMA  does  not  require  all  devices 
to  be  time-synchronized,  it  is  assumed  here,  for  simplicity  in  analysis, 
that  the  channel  time  axis  is  slotted  with  the  slot  size  equal  to  t  s., 
the  maximum  propagation  delay  between  all  pairs  of  users,  and  that 
all  users  are  synchronized  to  begin  transmission  only  at  slot  boun¬ 
daries.  The  CSMA  scheme  under  consideration  here  consists  of  the 

\ 

following.  A  user  with  a  packet  ready  for  transmission  (i.e.,  with 
a  packet  which  has  just  been  generated,  or  which  has  been  resched¬ 
uled  for  transmission  at  that  instant) ,  senses  the  channel  and  (i) 
if  the  channel  is  idle,  starts  transmitting  the  packet  at  the  be¬ 
ginning  of  the  next  slot,  and  (ii)  if  the  channel  is  busy,  then 
it  reschedules  the  transmission  of  the  packet  to  some  random  time 
in  the  future. 

We  consider  a  finite  population  of  M  users,  all  in  line 
of  sight  and  within  range  of  each  other,  such  that  each  user  can  be 
In  one  or  two  states:  backlogged  or  thinking.  In  the  thinking 
state,  a  user  generates  a  new  packet  (and  starts  transmitting  the 
packet,  if  the  channel  is  sensed  idle)  in  a  slot  with  probability 
o.  A  user  is  said  to  be  backlogged  if  it  has  a  packet  in  transmission 
or  awaiting  transmission.  It  remains  in  that  state  until  it  com¬ 
pletes  successful  transmission  of  the  packet,  at  which  time  it 
switches  to  the  thinking  state.  Thus  a  user  in  the  backlogged  state 
cannot  generate  a  new  packet  for  transmission.  The  rescheduling 
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delay  of  a  backlogged  packet  Is  assumed  to  be  geometrically  distri¬ 
buted,  i.e.,  each  backlogged  user  is  scheduled  to  resense  the 
channels  in  the  current  slot  with  probability  v;  as  specified  by 
the  protocol,  a  retransmission  would  result  only  if  the  channel  is 
sensed  idle.  In  this  model,  it  is  assumed  that  a  user  learns  about 
its  success  or  failure  instantaneously  at  the  end  of  its  transmis¬ 
sion  period.  All  packets  are  of  a  fixed  size  equal  to  T  slots. 
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III.  AVERAGE  STATIONARY  PERFORMANCE  [1,2] 

Let  n(t)  denote  the  number  of  backlogged  users  at  the 
beginning  of  slot  t.  We  observe  on  the  time  axis  an  alternate 
sequence  of  Idle  and  busy  periods  as  shown  In  Fig.  1.  We  follow 
the  approach  used  In  [1,2]  and  consider  the  imbedded  Markov  chain 
identified  by  the  first  slot  of  each  idle  period.  Using  properties 
of  regenerative  processes,  we  derive  the  average  channel  performance. 

Let  t^  and  t^r+^  be  two  consecutive  imbedded  slots; 
the  period  of  time  between  t^r^  and  t^r+^  is  called  a  cycle.  Let 
P  denote  the  transition  probability  matrix  between  t^r^  and  t^r+^ ; 
that  is,  the  (i.j)1"^  element  of  P^  is  defined  as 

?ij  =  Pr{n(t<r+1))  =  j  |  n(t£r))  =  i}.  0£i,j£M  (1) 

We  let  TP  denote  the  length  of  the  transmission  period.  If  the 

transmission  of  the  message  is  successful,  then  TP  =  T  +  1,  where 

the  additional  slot  accounts  for  the  propagation  delay  since  it 

is  only  one  slot  after  the  end  of  transmission  that  the  channel 

will  be  sensed  idle  by  all  users.  If  the  transmission  of  the 

message  is  unsuccessful,  then  TP  ■  T£  +  1,  where  Tc  £  T  is  the 

time  to  detect  the  collision  and  abort  all  transmissions  if  the 

collision  detection  feature  is  in  effect,  and  T  ■  T  otherwise. 

c 

(r) 

n(tfi  )  remains  invariant  over  the  entire  idle  period  I  (since 
according  to  the  CSMA  procedure,  a  new  arrival  sensing  the  channel 
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idle  would  transmit  with  probability  one).  See  Fig.  1.  Thus  for 

tC[t(r),  t^+I-l],  n(t)  »  n(t^).  Let  R  denote  the  transition 
6  6  6 

matrix  between  slot  t^+  I  -  1  and  t^  +  I.  Since  the  success 

e  e 

or  failure  of  the  transmission  is  a  function  of  the  number  of  users 

(r) 

becoming  ready  in  slot  t^  ,  we  write  R  as  R  *  S  +  F,  where  the 
(i,k)t*'  elements  of  £  and  are  defined  and  expressed  as 

sJt  =  Pr{n(t^+  I)  =  k  and  transmission  is  successful 
ik  e 

|  n(t(r)+I-l)  =  i) 
e 

0 

(l-a^tiuq-u)1'1] 
l-d-v)1  (l-a)M-i 

(M-i)CT(l-q)M~I~1(l-v)1 

i  M-i 

l-(l-v)  (l-o)” 

0  k>i+l 

(2) 


k<i 

k=i 

k=i+l 


7 


(l-g)M~i[l-(l-v)1-iv(l-v)i~1] 
l-(l-v)1  (l-a)M-i 

(M-i)a(l-a)M~1~1[l-(l-v) i] 
l-(l-v)1  (l-o)”'1 

M-i  \  ,,  ,M-k  k-i 

k-J  (1-°>  ° 

!-(!-») 1  (l-o)"'1 


k<i 

k=i 

k-i+1 

k-»i+l 

(3) 


During  the  transmission  period,  all  new  arrivals  join  the  backlog. 
Thus,  for  any  t€  [t^r^+I+l,  t^r^+I+TP],  we  let  jQ^denote  the  one- 
step  transition  matrix,  for  which  the  (i,k)tk  element  is  defined 
as  =  Pr{n(t)  =  k  |  n(t-l)  =  i}  and  expressed  as 

_j° 

qik  }(^)  (l-o)M_k 


k<i 

k>i 

w 


Finally,  to  represent  the  fact  that  a  successful  transmission 
decreases  the  backlog  by  1,  we  introduce  matrix  such  that 


its  (i,k)  element  is  given  by 

k-i-1 

otherwise  (5) 

The  transition  matrix  P  is  then  expressed  as 

T+l  T-+1 

P  *  S  5  j  +  f  £  c  (6) 

Let  n  »  {it.,  ir,  denote  the  stationary  probability  distri- 

bution  of  n(t^r^).  II  is  obtained  by  the  recursive  solution  of 

n  -  n  p. 

Since  n<t^  )  is  a  regenerative  process,  the  average 
stationary  channel  throughput  is  computed  as  the  ratio  of  the 
time  the  channel  is  carrying  successful  transmission  during  a  cycle 
averaged  over  all  cycles,  to  the  average  cycle  length.  Therefore 
we  have 


M 


E  *i  ».<«  t 

i-0 


M 

E’i 

i-0 


"i— r —  +  1  +  P  (i)T  +  (1-P  (i))Tc 
1-5 .  s  s 


(7) 
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i 

I 


where  P  (i)  is  the  probability  of  a  successful  transmission  during 
s 

a  cycle  with  n(t^r^)  ■  i,  and  is  given  by 


m  =  (M-^od-g)^1'1  (l-v)1  +  ivd-v)1"1  (l-o)**"1 

s  .  ..  .i  .M-i 

l-(l-v)  (l-o)  (8) 

and  where  (1-6^)  \  with  6^  =  (l-v)1  (l-o)M  is  the  average  idle 
( r) 

period  given  n(t^  )  =  i. 

Similarly,  the  average  channel  backlog  is  computed  as  the 
ratio  of  the  expected  sum  of  backlogs  over  all  slots  in  a  cycle 
(averaged  over  all  cycles),  to  the  average  cycle  length.  Therefore 
we  have 


n  = 


+  A(i) 


+  1  +  p  (i)T  +  (1-P  (i))Tc 
s  s 


(9) 


where  A(i)  is  the  expected  sum  of  backlogs  over  all  slots  in  the  busy 
period  with  n(t^r^)  =  i,  and  is  given  by* 


For  an  arbitrary  matrix  B,  we  adopt  the  notation  [B]  to 
th 

represent  the  (i,j)  element  of  B^. 
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T  M  Tc  M 

A(i)  “  EE  j[S  /ly  +  EE  j  [F  /] 

1-0  j-i  £-0  j-i 


M 


■  E  i 

j-i 


£=  0 


ij 


(10) 


By  Little's  result  [5],  the  average  packet  delay  (normalized  to  T) 
is  simply  expressed  as 

D  =  n/S  (n) 
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IV.  INTERDEPARTURE  TIME  DISTRIBUTION 


Theorem: 

The  generating  function  for  the  interdeparture  time  dis¬ 
tribution  is  given  by 

ID*(z)  “  A  [I  -  Pj  (z)]_1  P*(z)  H  (12) 

where 

i)  A  is  solution  of  A  =  A  “  Z  ^Tc+^]~^  A  S.T+^  J 


ii)  P^(z)  an<*  P^(z)  are  matrices  such  that  their  (i,j)t*1 
elements  are  defined  as 


[p*<A 

A  2T+1  Jij 

*  \  T+2 
(l-fi^z 

(13) 

ij 

1-6  jZ 

[s?(2>] 

[«vV 

(l-6i)zTc+2 

(14) 

ij 

1-6  ^z 

6±  =  (l-v)1 

\M-i 

(l-o) 

(15) 

iii)  H  is  a  column  vector  with  all  elements  equal  to  one. 
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Proof : 
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A  simple  recursive  procedure  exists  for  the  computation  of 


th  _(m)  A  dmID*(z) 

the  m  moment  of  ID.  Let  ID  -  - ^—L 

.  m 

dz 


,  ID.<»>U)  i  , 

- 


Z«1 


dz 


d“p*(z) 


pjC)U)  5  _-l_  4<»>W 

dz 


(a)/.  4  ^ 


dVkz) 


.  (We  use  the  convention 


dz 


that  the  derivative  of  a  vector  or  matrix  is  the  vector  or  matrix 
whose  elements  are  derivatives  of  the  corresponding  elements  in  the 
original  vector  or  matrix).  Clearly  E[ID]  =  ID^ ;  Var[ID]  *  ID^ 
+  ID(1)  -  [ID(1)]2,  ID(3)  *  E[ID(1D-1) (ID-2) ] ,  etc. 


Corollary: 

nj(m)  g^ven 

ID(m)  =  A  •  ID*(m)(l) 


where  ID*^ (1)  is  recursively  determined  by 


ID*(m)(l)  -  (I  -  Pj)-1 


m-1 


P*(m) (1) *H  +  (1) *H 


+  E  (")  pf(k?(D-i£*(m'k)a) 


(m-k) , 


k=l 


(21) 


(22) 


Proof : 

It  can  be  easily  proven  by  induction  that  differentiating 
Eq.  (18)  m  times  leads  to  the  following  relation 
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m*' 


ID*(m)(z)  -  P*(m)(z)H  +  (“)  p|(k)(z)*ID*(m~k)(z) 


,(m-k) , 


Letting  z=l  and  observing  that  ID*(1)  =  H,  and  P^-(l)  =  P^-  ,  we  get 


Eq.  (22). 


I.E.D. 


The  average  interdeparture  time,  in  particular,  is  given 


ID'"1'  =  A  (I  -  P^-)"1  j^P*(1)(l)  +  pj(1)(l)j  H 


±tl 


where  J^P^^(l)  +  P^^(l)^j  H  is  a  column  vector  whose  itk  element 


is  simply 


j=i-l 


iir-  +  T,  fit'  <T+1)  +  Pif  <Tc+1> 


The  variance  of  ID  is  given  by 


Var [ID]  =  A  (I-P^)-1  [p*(2)(1)  +  p|(2)(1)J  H  +  2P*(1) (1) ID*(1) (1) 


+  ID(1)  -  [ID(1)]2 


V.  DISTRIBUTION  OF  PACKET  DELAY 

Consider  an  imbedded  slot  t^r^  such  that  n(t^)  *  i, 

e  e 

and  let  D*(z)  denote  the  z-transform  of  the  distribution  of  delay 
(counted  starting  from  t^r^)  until  a  tagged  user  in  the  backlog  of 
size  i  is  successful.  Let  D*(z)  be  the  column  vector 
(D*(z),...,  D*(z))T. 

Theorem: 

D*(z)  is  given  by 

D*(z)  *  [I  -  Az)]"1  P*(z)  H  (26) 

—  —  — - s  — s  — 


where  H  is  the  column  vector  with  all  elements  equal  to  one,  and 
^  til 

P*(z)  and  P—  (z)  are  matrices  with  (i,i)  elements  defined  as 

“8  — S 


^(z)]ij 


(l-d^z 

l-5^z 


T+2 


(27) 


a™  i] 


ij 


(l-«i)z 

l-6^z 


T+2 


[*aTc+l]u 


.  Tc+2 
(l-6i)z  c 

1-6  * 


(28) 


[S  ] 


ik 


k-i 

otherwise 


(29) 
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(30) 


Proof : 

The  proof  is  similar  to  that  given  in  Section  III  for 
ID*(z) .  Noting  that 

[S  =  Pr{n(t^r^+I+1)  =  k  and  tagged  user  successful  I  n(t^r^)=i) 

S  XK.  G  6 

we  have 


"  T+1  (l-S  )zT+2 

D?(*)  ■  £  Ils  2  i'll  — ra 

j=i-l 


■5iz 


M 

E 

j=i-l 


,,  r  v  T+2 
T+1  (1_6l)z 

ij  — i rrr“ 


D*(z 


(31) 


hence  Eq.  (26). 


Q .  E .  D . 


Theorem: 

The  distribution  of  delay  is  given  by 


D*(z)  =  CqT+1)  zT+1  +  2  i(£>  I*(z) 

e-o 


(32) 


where  £qT+^  and  are  defined  in  Eqs.  (49-53)  in  the  proof. 


18 


ww .» snnpp  i  ?  f- ip 


Proof : 

To  complete  the  delay  calculation,  we  need  to  compute 

(t) 

,  the  probability  that  an  arbitrary  new  packet  arrives  in  a 

slot  which  is  t  slots  away  from  the  next  imbedded  slot,  and  finds 

Itself  at  the  beginning  of  this  Imbedded  slot  i.i  a  backlog  of  size 

k;  indeed,  its  delay  is  then  z^D*(z) . 

Let  t  and  t'  denote  two  successive  imbedded  slots 
e  e 

t'  *  t  +I+TP+1.  Let  a(t  )  be  an  indicator  such  that 
e  e  e 


if  no  new  arrivals  to  the  backlog  occurred 
in  the  cycle  [te>  t^~l] 

if  at  least  one  new  arrival  to  the  backlog 

occurred  in  the  cycle  [t  ,  t'-l]  (33) 

6  £ 


We  denote  by  t  an  imbedded  slot  such  that  a(t  )  =  1.  We  first  seek 
a  a 

the  stationary  distribution  of  n(t  ),  that  is  the  distribution  of 

a 

backlog  at  the  beginning  of  a  cycle  in  which  arrivals  occurred. 

Let  t^  and  denote  two  successive  imbedded  slots 

a  a 

such  that  a(t^)  =  a(t^r+^)  =  1.  The  imbedded  process 
a  a 

(n(t  ),  t  **...,  t^r\  t^r+^,...}  is  a  Markov  chain  with  transition 

probabilities  Pr{n(t^r+^)»  j  |  n(t^)=  i)  derived  as  follows. 

a  a 

Let  S- ,  F-,  Qr,  P—  and  P  be  transition  probability  matrices  defined 
0  0  *"*0  0  0 

00 


i 


1 

a 


I 

1 


1 
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l4]ij  “  Pr{n<c^  “  J»  “<**>  *  1  I  n<te)  “  i}  (40) 


(l-o)(M~i)(T+2)  fivd-v)1-1] 

i  M-i 

l-(l-\>)  (l-o)" 


J-i-1 


(1-g) (M~1) (Tc+2) [l-(l-v)1-lv(l-v)1~1] 


i  M-i 

l-(l-v)1  (l-o;  x 


otherwise 

(41) 


y±  '  Pr(a(te)  =  0  |  n(t£)  =  i)  -  fP^L  ^  +  t P^l ±  ± 


(42) 


Pr{a(t  )  =  1  |  n(t  )  =  i)  *  1-y. 
e  6  l 


(43) 


Let  P,,,  Piri,  I»  ,  and  P..  be  transition  probability  matrices  defined  as: 
— 11  — 10  — 00  01 


For  i  =  0,1,...,M-1;  j  =  0,1, 


[£ll]ij  “  =  j’  a(te}  =  1  1  n(te}  =  1’  a(te)  =  1} 


[P] 


ij 


(44) 
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•■».-'  c' itf  il;-.  "*~  '  ~  ■•ifcV'it. * 


[—10 5 i j  "  Prtn^P  =  “(tp  =  0  I  n^te>  ■  ®ue>  =  1} 


l^aJij  1-y, 


For  i  =  0,1, ...,M;  j  =  0,1,. . . ,M 


[Pooly  -  Pr{n(t;>  =  j,  a(t;>  =  0  |  n(te>  -  i,  o(te)  =  0} 


=  [Prl,, 


[^01]lj  "  =  3,  a(t;>  =  1  |  n(te)  =  i,  a(te)  =  0} 


It  is  clear  that  for  i,j  ^  M  we  have 


Pr{n(t<r+1))  =  j  !  n(t<r))  =  1} 


[-ll]ij  +  [-10(I  +  ^00  +  ^00  +”^  ^01]ij 


[-ll]ij  +  [-10  (I  “  ^01]ij 


Let  A  denote  the  above  transition  matrix  restricted  to  rows 
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i  =  0,1,..., M-l  and  columns  j  =  0,1,..., M-l.  Let  a  =  {ag,...,a^ 

(r) 

denote  the  stationary  distribution  of  n(t^  ).  £  is  obtained  as 

the  solution  of  A. 

Consider  an  imbedded  slot  t^r^  such  that  n(t^r^)  =  i  <  M 

a  a 

and  a(t^r^)  =  1.  Let  t^  denote  the  imbedded  slot  following  t^1^  . 

(£) 

Let  5^  (i)  denote  the  probability  that  a  tagged  packet  generated 

in  the  cycle  [t^r^,t'~l]  arrives  in  slot  t'  -  £  -  1  and  finds  itself 

at  the  beginning  of  t^  in  a  backlog  of  size  k.  We  use  the  value 

"k  =  0"  to  represent  the  case  where  the  arriving  packet  starts 

successful  transmission  at  its  arrival;  clearly  in  this  case  the 

(r) 

new  packet  must  have  arrived  in  slot  t^  +1,  and  £  must  equal 

T  +  1.  We  now  proceed  with  the  derivation  of  (i) . 

Let  m^,  0  <  <  M-i,  denote  the  number  of  arrivals  in 

slot  t*  -  £  -  1.  If  the  transmission  is  successful,  then  0  <  £  <  T+l 
e  —  — 

we  define  in  this  case  p  to  be 

v 


T+l 


j=T+l-v 


0  <  v  <  T+l 


If  the  transmission  is  unsuccessful,  then  0  £  £  <_  T£+l;  we  define 

in  this  case  p '  to  be 
v 


T  +1 


c 


j=T  +l-v 
J  c 


m. 

J 


0  <  v  <  T  +1 
—  —  c 
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For  pT+1  1  M-i  and  p'  +1  <_  M-i,  we  have 


Pl"{  *  *  *  ,mT+l)  *  transn,ission  successful 


n<tar))  =  i,  “(tfr))  =  1} 

u  a 


n 


A  -A  1+U  >i+P  .1 

0  v=0 _ y  y-f  1 

i-Yi 


(49) 


Pf(  (m0»mi>  •  •  ♦  +l^’  transmission  unsuccessful 


I  n(tgr))  =  i,  «(tfr))  =  1} 


(r), 

a 


111 i*i+U0  v  0  lq)i+K'i+Vv+l 


(49a) 


It  is  also  clear  that  for  0  ^  ^  ^  T  +  1  and  0  <  k  <  M,  we  have 


pr{ tagged  user  arrives  in  slot  t'-£-l  and  finds  itself 

e 

in  a  backlog  of  size  k  |  (ra^ra^ . . .  ) , 

transmission  successful,  n(t^r^)=i,  a(t^r^)=l} 
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—  n  *-m.  (k-o).  .T+1-i,  i<Ul+li«-i 


k-'i+l  if  °-^-T  *  i-k-M*  uT+l=k_i‘fl 


0  otherwise 


Similarly,  for  0  £  t  £  Tc+1  an<*  ^  —  k  —  we  kave 


Pr( tagged  user  arrives  in  slot  t^--£-l  and  finds  itself 


in  a  backlog  of  size  k  |  (m^.m^, . . .  ,mT 
transmission  unsuccessful,  n(t^r^)=i,  a(t^r^)=l) 


if  0<^<TC+1,  i+l<k<M,  \i'  -=k-i 

k.-l  fpil 


0  otherwise 


(50a) 


Therefore  we  can  write 


«ru«>-  z 


,  *-*i,i+l  nn  f^i+p  ,i+p 

1  v=0  v  \ 


<m0,  •  •  •  T+l 

PT+1  -  M_1 


•  E 


VT+1=1 


and  for  i  £  k  _<  M 

r 


4£>(i)  =  ■< 


ul 


(Mq>  •  •  • 
HT+1=k-i+l 


k-i+1 


ISll-1+»0  "o  I51‘\'I*W1 


1-Yi 

if  TC+2<£<T  and  i<k<M 


and  if  0<£.<TC+1  and  k=i 
T 

ml  “Ji,i+W0  v=0  l£L]i+\>i+\+l 


ts]. 


(in^  j  •  •  ♦ 

UT+l=k-i+1 


♦  E 

(iHq  >•••,111, 


k-i+1 


1~Yj 


"i  iVW0  "o 


k-i 


Tr+1) 


“ic+i  ‘  k‘I 


i-y< 


if  0<£<Tc+1  and  i+l<k<M 


which  reduces  to 
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W  «£  k*±?l  *-  S^'i.x  ‘S'x.x+m,  tAl^>k 

2~i  t_4+i  2-j 


.  k-i+l  . 
m^=l  x-i 


i-Yi 

if  TC+2<£<T  and  i<k<M 
and  if  0<£<TC+1  and  k=i 


k-i+l  _  k+l-m„  [S  ^T_'e].  [Q] 


l/j]. 


®£  Wt  * — *  'i,x  ^  x,x+tn^  -1  —  x+m^jk 

2-*  k-i+l  ^ 


.  k-i+l  . 
m^=l  x=i 


k-i 


1-Y, 


Tr+1-, 


c(*>(i)=!J  Cl?1»t4ini  t5  C  }i*“e»k 

1  .  _ i  k-i 


“r1 


l-Y, 


if  £=TC+1  and  i+l<k<M 


kyt1  k+^p£  ^  5i,x  ^x,x+n£  -^x+m^.k 

JLj  wTTTT  £—* 


,  k-i+l 

ttl£=l  X=1 


1-Y, 


k-1 


m^l 


k-i 


kg£  naT^i1>lt  i2ix>lrt.t 

x=i 


^i 

if  0</<Tc  and  i+l<k<M 
(52) 


We  remove  the  condition  on  n(t^r^)  by  noting  that  =  Pr{n(t^r^)=i 

j  a(t^r^)=l),  and  thus  write 
a 


26a 


(53) 


E 


i=0 


Let  5.^  ^  -  (£^\  » • •  • »  The  distribution  of  packet 

T+l  (T+l}  P 

delay  is  z  with  probability  £q  ,  and  z  D*(z)  with  probability 

(Jf) 

5^ hence  Eq.  (32). 

Q.E.D. 

As  with  inter departure  times,  a  simple  recursive  procedure 
exists  for  the  computation  of  all  moments  of  the  delay. 


Let  D(">  5  CS ii  ,  „*<"><*>  5  -'•><»)  P. <»)(,)  S 

j  m  *  -  .  m  *  v  7 


dV^z) 


dz 


z*=l 


dz 


dz 


m 


and  P*(m)(z)  = 


dV^z) 


dz 


m 


Corollary; 

is  given  by 


(m)  _  (T+l) !  (T+l) 

(T+l-m) !  C0  + 


T 

E 

1*0 


mm 

(£,m) 


E(;)^iltv-k)«. 

j_k=0 

(54) 


where  D*^m^(l)  is  recursively  determined  by 
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/ 


D*(m)(l)  -  [I  -  P-(l)]"1 
'  8 


P*^(1)H  +  p4(m)(l)H 

— S  —  —8  ~ 


m-1 

E  (*)  pj(k)a)B*(”"k)a> 

k=l 


(55) 


Proof: 

By  successive  differentiations  of  Eq.  (31)  and  Eq.  (32), 
and  letting  z=l,  we  can  easily  establish  Eqs.  (54,55). 

Q.E.D. 
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VI.  CONCLUSION 


We  derived  in  this  paper  the  distributions  of  packet  delay 
and  interdeparture  time  for  CSMA  channels  with  finite  populations  of 
interactive  users.  The  approach  is  similar  in  nature  to  that  used 
for  slotted  ALOHA  in  [3],  except  that  the  basic  model  underlying  the 
analysis  is  here  semi-markovian  rather  than  markovian.  We  have 
derived  simple  recursive  procedures  to  compute  all  moments  of  packet 
delay  and  interdeparture  time,  and  closed  form  expressions  for  their 
generating  functions. 

Contrary  to  their  counterpart  in  slotted  ALOHA,  the  gener¬ 
ating  functions  derived  for  CSMA  may  not  prove  very  practical  in  the 
numerical  computation  of  the  distributions.  Indeed,  they  require 
symbolic  inversion  of  matrices  whose  elements  are  themselves  z-trans- 
forms.  However,  a  numerical  procedure  which  allows  to  compute  approx¬ 
imations  to  the  distributions  can  be  devised  as  follows.  Consider  first 
the  case  of  interdeparture  time.  It  is  easily  realized  that 


ID*(z) 


00 


m=l 


where 


'  (ID8,<«ow 


,IDM,  (m) 


(56) 


(57) 


I°?m)  “  ^f(z)  — tm-1)  (z^  for  m  i  2 


(58) 
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Eqs .  (58,59)  are  equivalent  to 


M 


”>?,(  »<*> 


-  _S[id] 


(60) 


j=i-l 


and  for  m  >  2 


M 


ID* 


,  ,(z)  -  V  Pj  (1-6.)zTc+2  1+S,z+62z2+...  ID* 

« (®)  (_-dJii  i  i  1  J  j. 


J-i-1 


(m-1) 

(61) 


Thus  by  successive  polynomial  multiplications  and  additions  one  can 
generate  numerically  an  approximation  of  the  distribution  of  ID,  the 
accuracy  of  which  is  a  function  of  the  position  at  which  the  infinite 
series  are  truncated  and  of  the  maximum  value  given  to  m.  A  similar 
procedure  can  be  devised  for  the  distribution  of  delay. 
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